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1.0 INTRODUCTI ON ^ ^UALTfY 

Th« gravitational gaopotantial la uaually rapraaantad by tha apharioal hamonic 
axpanaion solutior to Laplaoa'a aquation. Tha haraonlo ooaffioianta ray ba da- 
tamlnad by raduotion of aatalUta trajaotory data, naan ooaan aurfaoa altlaa- 
tar raaaurananta, or local aQcalar<w«atar raadinga 2), Tha aonal harronics 
giva rlaa to aaoular variations in tha angular orbital alaaanta but only to pari- 
odio ohangaa in tha roiaanta alarants. Tha sacular ^nd aactorial harmonics 
oausa no trua aaoular variations, but tha non-taro naan of the pariodic varia- 
tions in the raar motion Irduoas a linear perturbation in tha mean anomaly (3). 

Tha nomallxad hamonio ooaffioiants do not daoraasa in magnltuda as the degree 
of the expansion increases. Also, each harmo.tlc is praaultlplied by a pover of 
tha ratio of tha Farth radius and tha position radius. Since for low Earth sat- 
ellites this ratio 1s near unity, successive terra in tha axparsicn do not tend 
to diminish. This necessarily requires the inclusion of a vary large number of 
terra to accurately model the gravitational perturbations. Recursion elgorlthrs 
give an efficient means for calculating tha higher degree perturbations, but 
the computation times grow geometrically with the degree (H, The compu- 

tational environment may preclude tha use of such expensive models, and the 
common trade-off is to simply truncate to only a few Important terms. 

The ooafflolents of the large degree models are determined from observations as 
a complete set, and it Is ret clear that simply truncating the model Is the 
most appropriate action. Certainly, fitting the observations with only the 
truncated model v'ould yield different numerical values for the harmonies. The 
intent of this study Is to develop a fitted truncated model and analyze any dif- 
ferences between tMs **fltted" model and one derived by simply trunontlrg. 

Based on the study, recommendations are made for an appropriate model for use 
In a mission planning environment. 
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2.0 ERROR NORMS 

Pitting Qod«ls iitapli«8 a nors by triiioh to maaaura the quality of the fit. The 
harmonioa are gener;'ted by minimizing the residuals bet«reen the observations 
and the model predictions. Observations may be very Indirect measurements 
of the model and may Include errors and biases so that minimizing the resi- 
duals becomes no easy chore. To simplify natters a bit, assume that a 
reference geopotential model represents the real world exactly. This truth 
model then can be used to evaluate the accuracy of any other model and eliminate 
many difficulties resulting from real world observations. The most appropriate 
norm should measure the differences In the state vectors defined by the equa- 
tions of motion derived from the different gravity models. Suppose also 
that an^exact solution of the equations of notion Is available. If one de- 
fines x(xo,Vo*^^ position vector at all times t for the truth 

model such that 


X<Xo»Xo* 


t=0) = 


( 2 . 1 ) 


and a similar solution to the truncated model, then the norm 



*<Xo.Xo.t) 


f 


dtdxodxo 


( 2 . 2 ) 


gives a measure of the differences between potential models. Here, the 11*11 
is the dot product norm and the domains fl and represent all possible 
values of the initial state vectors x^ and v^. Such a norm is a good 
indicator of the propagation errors resulting from different models, but 
it is not an easily computable norm. More direct norms on the geopotential 
model itself may be a worthy substitute. For Instancj;, the norm of 

the difference between the potential models V and V is: 

Jr [ (V - To* dx' ■ (? 3 ) 


For this norm^to be a minimum, the first variation must vanish so that the 
error e = V-V roust satisfy 

I eu)dx s 0 V U) e L* 


2 


( 2.^0 
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The test function w nay be thought of as the first variation of the potential 
V, and represents the space of all functions which are square Int^rable. 
Only a finite nunber of the basis functions are available in the truncated 
model so that the error can be tested to only a finite number of functions 
0 ). 

Since it is the gradients of the potential' and not simply the potential 
which define the accelerations, a more appropriate norm may be the semi- 
norm 



(2.5) 


In this norm, the error must be orthogonal to every test function in the 
space of functions tdiose derivatives are square integrable. 


I • 


V(D dx = 0 V u) e 


fi 


( 2 . 6 ) 


In spherical coordinates r, y, X the above equation reads 


3o) ^ (u^ - 1) ^ ^ 1 ^ ^ 

hr 3r r^ 3y 3y r*(1 - y*) 3X 3X 


drdydX = 0 


(2.7) 


This formulation is singular at the poles, but this may be easily averted 
by multiplying through by r*(1-y*) to give 



2/4 2 \ 3c /4 2 \^ 3® 3u) ♦ 3e 3u) 

'*(1 - y *) — — ♦ (1 - y*) — — rr — 

3r 3r 


3y 3^ 3X 3X 


drdydX s 0 


( 2 . 8 ) 


which corresponds to the minimizing of the weighted norm 


J 


= i r^(^ - P*)Ve . Ve dx 


n 


(2.Q) 
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3.0 THE GEOPOTEHTIAL EXPANSION AHD TF5T FUNCTIOHS 

A 

The function V ant V are represented by the usual spherical hamonio 
expansion 


V s - r (R/r) Z P]>] ||)(vi){C|>)]DCOSinX ♦ SfjiuSinDX} 
r nsO m=0 ’ 


(3.1) 


and 


. N n n 

^ ' r Lo SnjjSinmX} 


(3.2) 


trtiere 


Cnm * (^nm 


(3.3) 


and 


^nm ~ ^nm 


^Sn„/P 


nm'*^nni 


(3.10 


Here P^^ Is the average of the Legendre functi e over the entire latitude 

A A 

domain. C|]q] and are given by the truth model, and the unknowns 

in the truncated model are the normalized harmonic deviations. Solving 
for the normalized deviations, Instead of simply the harmonic coefficients 
themselves, conditions the matrices of the linear algebraic problem resulting 
from either equation (2.^) or (2.6). 

The test functions are then 


0] = 



^ ,iu 

Pki 


k = 0,N 
I s 0,k 


(3.5) 


and Its associated derivatives are 


= -(k+1) /R\k 
8r r^ \r/ 


Pm 


(3.6) 
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l .iix 

Sy r \r/ dy / 


and 


(?.7) 


3b) 

M tm 

3X 


r \r/ Pyp 


( ■ 


In spherical coordinates, the iBinimuin of the L* norm must satisfy 


.R+h 


I I I e b) dr dX dy s 0 

-e \ ■'b 


(3.9) 


where h defines the radius domain and e the latitude domain. 


iAj 
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H.O THE LIWEAH DISCRETE EQUATION?! 

Replacing the error and teat functions Into equation (?.4), one arrives at 


c e R+h ( N 


L .! i jto 
.1 


y y 

msO 


1 /"X"*" VKt 
Wkl 




ISSninSinmX 


e^^l drdpdX 


♦ SnnSinmX 


,lix( 


drdydX 


ZCnigCoaaX 


c a R+h ( N n 1 /*’V 

/ i / S S ;« (;) 

-C R ( n«N-«-1 a*0 


PnaPkl 


PnmPkl 


Cnncosmx 




but, duo to the orthogonality of the transcendental functions, these equations 
ray be oonsi<'erably reduced. Defining 


R+h 1 / R\n+k 

/ ;«(;) 


dr 


(14.?) 


pn,m _ 


/ Pnm Pkm 


(4.3) 


and 


s D 
"k 


I kV 


^Pkn Pmn^ 


(4.4) 
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th« linear equations read 


N 

T. 

nsm 




nm 


A 

N 

I 

nsN-Kl 




(H.5) 


N 

Z 

nsm 




nm 


A 

N 

£ 

nsH+1 


*n,m 

*k 




(H.6) 


for m s 0, . . . N 
k s m, N 


This is A system of linear algebraic equations In the unknowns ^nm 
®n«- Only the harmonics of the same order are coupled.^ gor the ease where 
the fit is applied over the w^lOle domain c s l, then a"* s 0 except 
when k 8 n. In such a case, the problem is completely uncoupled and the 
harmonic deviations must vanish. It Indicates that, in this particular 
norm, the simply truncated model io Indeed the optimal fit. This is not 
the case for restricted domains or for other norms, as will be shown. 


The weighted semi-norm also results in a linear sysvem of equations. 

With the definitions 


pn,m 

k 




'^^km 

dy 


<^Pnm 

dy 


dy 


(4.7) 




W^)Pkm ^nm 


(4.P) 


and 


n,ra 

'‘k 


(k ♦ 1)(n ♦ Dg"’® + F 


,r.m 




(PkH^nm) 


the equations still retain the same form as those in equations (4.5) and 
(4.6). Since the derivatives oi the Legendre functions are not orthogonal, 
these equations r«nain coupled even when fitting over the entire domain. 
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Th« linear system of equations may be solved using a Cholesky decomposition 
since they el's positive-definite and ayismetric. 
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5.0 THE QUADRATURES 

The quadraturea involving the radius aay be ooi^uted analytioally. Defining 
the ratio r, « 1/(1 * h/E), then the parameter is given by 


dJ = 1 (1 - (5.1) 

R(1 ♦ n + k) 


The quadraturea involving the Legendre functions may not be found analytioally 
for e ^ 1. A numerical Gauaf’ point rule is a convenient and arbitrarily accu- 
rate method to evaluate the quadrature. An s^^ order rulb is given by 


b 

f b - a 

I f(u)du = — 

a 


8 

I (ij®f(y?) ♦ 
i=l ^ ^ * 


(5.2) 


where 


a 


V 


i 


b - a b -f a 



(5.3) 


§® = GausL points 
til® = Gauss weights 


The residual is given by 


(b-e)^^®^^>(S!)'‘ 

R s 

® (2s ♦ 1)(?al)-^ 


f^^®^(y) ; ue(a,b) 


(5.’0 


Even though only s function evaluations are needed, the rule is accurate 
to order 2a -f 1. The dotuain defined by c nay be partitioned into smaller 
subdomains to take advantage of the factor 


q 
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(V)"" 

Recursion foniulas may be used to evaluate the Legendre polynomials and their 
derivatives at the Gauss points. These include 


(n - m ♦ 1)PrH.1,m * Pnm * “)Pn-1,m 


(?.5) 


and 


dP 


nm 


(Ii2 _ 1 ) r np P„B - (n ♦ m)P„.l n 

dp 


(5.6) 


with starting values 


Pam = O “ • 3 • 5 ... (2m - D) 


am 


(5.7) 


and 


^nKf-1,m ~ P mm 


(5.8) 


Lastly, the mean values of the Legendre polyncMnials are 


no - 


1 


/2n + 1 


(5.9) 


and 


nm 


(n * m) I 

?(n - m)!(2n + 1) 


m ^ 0 


( 5 . 10 ) 
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6.0 NUMERICJU EVALUATION OP THE FITTED TRUNCATED MODELS 

A program has baan davalooad to fit tba trunoatad modal to a rafaranoa modal, 
in this oasa tha GEM10 [ll. Fit is aooording to tha L^ or tha waightad h' 
semi-nom or a oombination of tha two. A ninth ordar Gauss intagrstion nils 
on two subdomains is amployad. To avaluata tha quality of tha fit, anothar pro- 
gram has been developed to propagate satallita orbits in a gravitational field. 
The recursion formulation desorlbad in Mueller (4] is Miployad to oomputa tha 
accelerations; and an embedded Runga-Kutta method (6) is used to numarloally 
integrate the cartesian eouatlons of motion. 

The state vector resulting from the fitted truncated model is compared to the 
state resulting from th^ truth model. As a control, the simply truncated model 
is also evaluated in this manner. The L^(T) error norm has been chosen to in- 
dicate the performance 


T 

j Vo. t) - x(xo, Vf„ t)||^dt 

0 


A variety of initial conditions can be used to test the models thoroughly. The 
time interval has been selected as one day for all test cases. 


6.1 TEST CASE 1 

In the first ease, a i»x4 model is fit to a reference 8x8 model. The fit la 
over the entire latitude domain and is up to 500 km in height. The norm used 
is the linear combination of the L^ and weighted semi-norm. The numeri- 
cal values of the fitted and simple normalized harmonics are shown in table 6-1. 

The initial conditions are such that the satellite remains in the 5C0 km height 
band. The height of the semi -major axis is initially 300 km. The eccentricity 
and argument of perigee have the nominal values of e s 0.01 and u s 0. The 
inclination and ascending nodes have been varied to scan a number of initial 
conditions. The inclinations range from 0 to 30 degrees in increments of 6 
degrees, while the node ranges from 0 to 360 degrees in 45 degree increments. 

The error norms for both the fitted and simply truncated models are shown in 
table 6-11. The results indicate that the fitted model shows little, xf any, 
improvement over simple truncation. 
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TABLF 6-1.- FITTED AND SIMPLY TFUNCATED NORNALI7ED HARMONICS (xIO*^) 


0 

1 -4PU.0791 

O.oi?6 


1 -484.1654 

I 

O.Q504 

1 

1 

I -0.0312 

2.0501 


1 0.0010 

2.0206 

2 

1 

1 2.4456 

O.59«^0 


1 2.4340 

1 

0.0027 

3 

1 

1 

0.7950 


I 

0.7003 


K 



1 

B 

1 

1 2 
1 


T 


1 

0.6604 

! 


! 

0.5411 

I 


I 


1 


1 

-0.5011 

1 

1 

1 -0.0542 

-0.5352 

1 


! -0.0024 


1 


1 

0.3313 

! 

2 

1 -1.3^60 

0.3521 

1 


1 -1.3991 


1 


1 

1.0452 

! 

3 

I 

O.Q60e: 

1 


! 


1 


1 

-0.1701 

1 

4 

1 

-0.1953 

1 


1 


1 


! 


I 

3 4 I 

I 

I 

I 

I 

0.3061 -0.4957 I 

0.2520 -0.460? I 

I 

-0.5279 0.P?27 I 

-0.6235 0.6640 I 

! 

1.4546 -0.1976 ! 

1.4125 -0.2018 I 

I 

0.3734 I 
0.2908 I 
! 


GM « 3.98600657x101** (m? s“2) 
3.90600640x101** 


TAPLE 6-II.- POSITION ERROR NORMS (FITTED/SIMPLE) 


Inclination (deg) 



1 

1 

ms' 

0 1 

6 1 

12 1 

10 1 

24 1 

30 


i' 

1 

1 

0 

1 

3. 6/4.0 ! 

1 

3. 3/4. 6 ! 

2. 0/4.0 1 

2.3/3. 3 J 

2.0/2. 8 1 

1.9/2. 6 

n 

1 

1 

1 

45 

3.‘=/2.6 I 

3.0/2. 3 I 

2. 4/1. 6 ! 

I 

1. 8/1.0 I 

1. 5/0.0 ! 

1. 5/1.0 

0 

d 

1 

1 

I 

QO 

4.7/4. 8 I 

5. 3/*'. 5 1 

6. 2/6.6 1 

6. 9/7. 6 1 

7.2/8. 1 1 

6.9/8. 1 

% 

1 

1 

t 

135 

6. 4/8.0 I 

5.<>/7.4 I 

5.0/6. 3 ! 

4. 0/5.0 ! 

3.5/4. 2 1 

3. 6/4.0 

G 

e 

1 

1 

t 

100 

1.9/1. 6 ! 

2. 0/1. 7 I 

2.3/2. 1 1 

2.0/?. 8 I 

3. 3/3. 4 I 

3. 6/3. 9 


1 

I 

1 

225 

5. 1/5.8 ! 

5. 1/5.6 I 

4.8/5.? 1 

4. 3/4. 8 ! 

3.7/4. 1 ! 

3. 2/3. 4 


1 

I 

1 

270 

6. 2/7. 3 ! 

6.0/7. 1 ! 

5. 3/6. 3 ! 

4.2/5. 1 ! 

3.5/4. 2 1 

» 1 

3. 3/3.8 


1 

1 

I 

31*^ 

6. 2/5.7 1 

6. 4/6.0 1 

7.3/7. 1 1 

8. 6/8. 6 ! 

0. 6/9.0 ! 

O.O/1O.0 


1 ? 
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6.2 TEST CASE 2 


ORiaWAL PAG£ B 
Qp O 


The second test case again ases a aodel fit to an 8x6 referenoe. This tiM 
the latitude domain is restricted to 30 degrees for evaluating the tesseral and 
sectorial harmonic fit. The zonal harmonics, however, use the full 90-degree do- 
main. The results of the fit are shown in table 6-III and the position error 
norms resulting from the fit are glv«i in table 6-IV. Again, the fitted model 
shows no Improvement over simple truncation. 


TABLE 6 -I II.- FITTED AND SIMPLY TRUNCATED NORMALIZED HARMONICS Xx10”6) 


1 

* ^nm 

2 

3 

4 

1 

1 

Smb 

2 

3 

4 

! 0 

-1184.0791 

0.9136 

0.6604 

T 

1 






-484.1654 

0.9584 

0.5411 

1 

1 

1 





1 1 

-0.0276 

2.5849 

-0.2356 

1 

-0.1072 

0.6420 

-0.4325 


0.0010 

2.0286 

-0.5352 

! 

! 

! 


-0.0024 

0.2520 

-0.4693 

! 2 

2.5043 

0.3371 

0.3866 

2 

-1.0655 

-0.2916 

1.2699 


2.4340 

0,8921 

0.3521 

! 

! 

1 


-1.3951 

-0.6235 

0.6640 

I 3 


0.9593 

1.0819 

3 


1.5002 

-0.2611 



0.7003 

0.9885 

! 

1 

! 



1.4125 

-0.2018 

! n 



-0.1588 

4 



0.4852 




- 0 . 1953 

1 

I 




0.2988 


GM = 3.98600657x10 (m3 s“2) 
3.986006»I0x10 
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TABLE 6-IV.- POSITION ERROR NORPS (FITTED/SIMPLE) 


ms I 


inollnatlcNn (dag) 


12 


18 


2«l 


30 


n 

o 

d 

e 

d 

e 

g 


! 

1 0 

1 

I 

1 

1 

3.7/tl.9 

1 

1 

• 

3.5/11.6 

I 

1 

1 

3.0/A.0 

I 

1 

1 

2. 4/3. 3 

r 

I 

1 

2. 0/2. 8 

1 

1 

1 

1.8/2. 6 

1 

1 H5 

I 

I 

a 

11.9/2.6 

1 

I 

1 

11.7/2.3 

1 

! 

a 

4. 0/1. 6 

1 

1 

1 

3. 3/1.0 

1 

1 

1 

2. 7/0.8 

1 

1 

a 

2. 4/1.0 

1 

1 90 

I 

1 

I 

1I.7/H.8 

1 

! 

1 

5. 1/5.5 

1 

! 

t 

5. 8/6. 6 

1 

1 

1 

6. 5/7. 6 

1 

! 

t 

6.7/8. 1 

I 

1 

I 

6.6/8.’ 

1 

1 135 

• 

1 

1 

5. 0/8.0 

1 

! 

1 

H. 7/7. 11 

1 

I 

1 

4. 0/6. 3 

1 

1 

1 

3. 2/5.0 

1 

I 

a 

2. 8/4. 2 

1 

! 

a 

3. 0/4.0 

1 

! 180 

1 

! 

1 

2. 0/1. 6 

1 

! 

1 

2. 1/1. 7 

1 

1 

1 

2.6/2. 1 

1 

! 

1 

3. 3/2. 8 

1 

! 

a 

3.9/3. 4 

1 

1 

t 

4. 3/3. 9 

1 

! 225 

I 

1 

1 

1 

5. 7/5. 8 

1 

! 

t 

5.5/5. 6 

1 

1 

1 

4. 9/5. 3 

1 

! 

1 

4.2/4. 8 

1 

1 

a 

3.4/4. 1 

1 

1 

a 

2. 8/3. 4 

1 

1 270 

1 

1 

I 

1 

5. 5/7. 3 

1 

! 

i 

5.2/7. 1 

1 

1 

1 

4. 4/6. 3 

1 

I 

f 

3.5/5. 1 

1 

I 

1 

2. 9/4. 2 

1 

! 

1 

3. 2/3. 8 

1 

! 315 

1 

! 

1 

! 

6. 5/5. 7 

1 

1 

! 

6. 7/6.0 

1 

1 

1 

7.5/7. 1 

1 

t 

1 

8. 5/8. 6 

1 

I 

f 

9. 2/9. 9 

1 

! 

1 

9.2/10.0 
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7.0 BFC0>«HRNDj>TI0N.*! AND COKCLUSIONS 


OF 


poo*^ 




From the nuRerlcal results, it appears that there is no advantage to fitting 
truncated isodels (at least in the norms Investigated). If trunoation is a ne- 
cessity in the computational environment, the question left open is at what de- 
gree truncition should occur. The magnitudes of each term In the expansion can 
be approximated by the normalized harmonies multiplied by the distance ration 
power or 


J J- > 


(7.1) 


At a height of 300 km, only 4 terms have magnitudes greater than 10“^. They are 
n=?, ir=0; ns2, ms2; ns3, msl; and n=3, ms3. If the toler- 
ance is increased by one digit to 10“7, then the number of additional terms in- 
creases dramatically to 33, or most of the terms of a 10x10 model. 

The analyst also should consider other factors in selecting a model. Since in 
low Earth orbits the downtrack quadratic drag errors overi^elra the linear down- 
track errors due to truncation, it is not important to include the tesseral and 
sectorial harmonics. These terms contribute secular variations only in the 
downtrack. The zonal harmonics, on the other hand, contribute radial and cross- 
track secular perturbations which are much larger than the corresponding linear 
errors due to drag. A Iw degree model appears to be a very suitable model in 
the presence of strong drag (h < 300 km). In these strong drag cases, reso- 
nance is not an Important factor since drag inhibits the period from resonating 
with any particular harmonic over any significant length of time. For the weak 
drag cases, the analyst should include only that particular resonating term. 

Put recursion algorithms for computing the accelerations must build up a table 
starting from a low degree, so selecting a particular term is not always possi- 
ble. 

Another consideration in selecting a model is consistency and uniformity. The 
analyst should select models wh’'.ch are consistent with the realtime environ- 
ment. For instance, if a high degree of accuracy is required, one should cer- 
tainly select the same model which was used to generate the initial state from 
the observations. The analyst also may wish to select a simple model which is 
consistent with those used by other analysts in related work. In this way, a 
source of discrepancies can be eliminated. 
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